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ABSTRACT 


An  experimental  and  theoretical  study  has  been  made  oil 
the  history  of  the  pollutants  carbon  monoxide  (CO) ,  unburned 
hydrocarbons  (HC)  and  nitrogen  oxides  (NO  )  in  a  turbojet 
afterburner.  Experimental  traverses  at  several  axial  stations 
were  performed  in  a  simulated  afterburner  in  which  exhaust  from 
a  J-47  combustor  can,  operated  at  medium  power,  was  mixed  with 
fuel  spray.  Experiments  were  carried  out  both  in  a  non-bypass 
and  in  a  bypass  configuration  (secondary  air  was  mixed  with 
primary  exhaust) .  The  non-bypass  tests  were  carried  out  at 
high  combustor  efficiency,  and  yielded  the  following:  CO  = 

300  ppm,  HC  less  than  10  ppm,  NO  -  100  ppm.  In  the  bypass 
tests,  fuel  distribution  was  nonfinlform  and  combustor  efficiency 
was  low.  The  concentrations  obtained  were  CO  =  10,000  ppm,  HC  = 
1000  ppm,  NO  »  100  ppm.  The  theoretical  analysis  consisted  of 
a  computer  program  for  reacting  flow  with  turbulent  mixing. 

The  computer  program  was  very  slow  and  therefore  of  limited 
usefulness  in  terms  of  cost  and  questionable  results,  since  it 
could  not  be  checked  against  experiment  Infrared  measurements  of 
NO  in  the  combustion  tunnel  were  attempted.  Indications  were 
obtained  of  NO  at  the  5.3  micron  band,  but  quantitative 
measurements  were  not  obtained. 
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SECTION  I 


INTRODUCTION 


This  users  manual  was  prepared  under  contract  number  F33615-71-C-1125, 
and  represents  Item  Number  0001  Sequence  Number  A002  under  require¬ 
ments  of  that  contract.  The  purpose  of  writing  the  computer  pro¬ 
gram  was  'to  provide  the  Air  Force  with  a  means  of  estimating  air 
pollution  concentrations  in  the  exhaust  gases  from  afterburning  tur¬ 
bojet  engines. 
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SECTION  II 


COMPUTER  PROGRAM  SYSTEM  CAPABILITIES 


2  . 1  PURPOSE 


The  purpose  of  this  computer  program  is  to  solve  the  partial 
differential  equations  of  change  for  momentum,  species,  and 
enthalpy,  to  obtain  radial  and  axial  profiles  of  velocity, 
temperature,  and  species  concentrations  in  a  turbulent  reacting 
fluid.  A  finite  difference  approximation  to  the  equations  is 
made. 

The  chemical  species  considered  and  chemical  reactions  involved 
are  those  pertinent  to  burning  of  a  liquid  hydrocarbon  fuel  in 
a  hot  exhaust  gas  from  a  primary  combustor.  By  controlling  the 
input  a  calculation  of  the  mixing  and  combustion  in  a  bypass  type 
engine  may  also  be  made.  While  the  primary  purpose  of  the  re¬ 
search  has  been  to  estimate  levels  of  the  pollutants  NO,  CO,  and 
unburned  hydrocarbons,  in  exhaust  gases,  the  levels  of  several 
other  species  must  be  taken  into  account  in  order  to  calculate 
pollutant  levels. 


2  .2  GENERAL  DESCRIPTION  OF  THE  SYSTEM 


A  FORTRAN  program  has  been  written  to  achieve  the  purpose  descri¬ 
bed  above.  The  differential  equations,  transformed  frcm  (length, 
radius),  to  (length,  stream  function)  coordinates,  have  been  appro¬ 
ximated  by  finite  difference  equations. 

Initially,  certain  thermodynamic  data  for  the  chemical  species  are 
read,  followed  by  data  necessary  to  calculate  chemical  reactijn 
rates.  Then,  initial  data  of  temperature,  velocity,  and  species 
concentrations  as  a  function  of  radial  position  in  an  axisymmetric 
flow  are  read  in.  The  program  does  not  read  any  more  data. 

From  the  initial  profiles  as  a  function  of  radius,  a  conversion  to 
the  stream  function  variable  is  made,  with  appropriate  interpolation. 

An  iteration  scheme  is  used  to  solve  for  dependent  variables 
(velocity,  enthalpy,  species  concentration)  as  a  function  of  axial 
distance.  The  program  proceeds  stepwise  from  the  initial  condition, 
until  a  maximum  axial  distance  is  reached.  The  outer  boundary  con¬ 
dition  is  an  inviscid  reacting  flow.  When  the  flow  variables  at  the 
stream  function  location  adjacent  to  the  inviscid  flow  deviate  from 
the  inviscid  flow  by  1  part  in  1C^,  the  viscous  flow  region  is  ex¬ 
panded  by  one  step  in  stream  function.  The  number  of  stream  function 
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steps  is  allowed  to  increase  to  50,  if  necessary.  If  more  expan¬ 
sion  is  needed,  the  stream  function  step  size  is  doubled,  and 
number  of  steps  is  halved.  This  process  repeats  itself  as  necessary. 


2.3  FUNCTIONS  PERFORMED 


There  are  13  subdivisions  of  the  program. 

1.  afAIN  Program 

2.  Subroutine  SPECIE  (KK1,  KK2 ,  NU) 

3.  Subroutine  THERMO  (KRIS,  JA,  IT) 

4.  Subroutine  INIT 

5.  Subroutine  TRANS P  (NLOW,  NHIGH) 

6.  Subroutine  INTEG  (NLOW,  NHIGH,  IN) 

7.  Subroutine  SETUP  (RUP) 

8.  Subroutine  SETYP  (RUP) 

9.  Subroutine  SETHP  (RUP) 

10.  Subroutine  GAUSS 

11.  Subroutine  CHANGE  (RUP) 

12.  Subroutine  PRINT  (KUSK,KUN) 

13.  Subroutine  CHEMIE  (KPYS,  KRAS,  IJ) 

The  MAIN  program  controls  the  flow  of  information,  causing  THERMO, 
CHEMIE,  and  INIT  to  read  in  initial  data.  It  then  proceeds  with 
the  computation  and  prints  out  results  at  every  50th  axial  station, 
until  stopped  by  either  reaching  the  maximum  desired  axial  distance 
(ZMAX  feet)  or  exceeding  the  time  limitation. 


2.4  FUNCTION  DESCRIPTIONS 


2.4.1  Main  Program 


Summary  of  Operations  performed  by  MAIN 


The  main  program  first  sets  certain  parameters,  namely  ISP  =  14, 

IRE  =  17,  IMU  =  11,  which  ara  in  labelled  common  ZZZ.  There  are, 
respectively,  the  number  of  species,  the  number  of  reactions,  and 
the  number  of  the  last  reaction  which  involves  three  bodies  (seeo 
CHEMIE).  In  addition,  the  variable  TAU  is  ^et  equal  to  .4  x  10  , 

which  is  the  hydrocarbon  burning  time  in  seconds.  TAU  is  in  labeled 
common  XXX.  Certain  other  parameters  are  also  set,  so  as  to  cause 
THERMO  to  read  thermodynamic  data,  and  CHEMIE  to  read  chemical  reac¬ 
tion  rate  data.  The  main  program  then  calls  upon  INIT  to  read  in 
velocity,  temperature,  and  species  concentrations  as  a  function  of 
radius.  The  first  item  actually  read  is  NNN,  which  defines  the 
number  of  radial  locations  considered.  INIT  converts  the  indepen- 


dent  variable  from  radius  to  stream  function,  and  interpolates  to 
obtain  initial  dependent  variables,  velocity,  temperature,  and 
species  concentrations,  at  equal  intervals  of  stream  function. 

The  stream  function  interval  is  denoted  by  DPSI.  DELX  is  the 
interval  in  feet,  and  is  defined  to  be  .00001.  RUP  =  DELX/DPSI**2 
is  defined.  XMAX  =  0.7  (feet)  is  defined  as  the  maximum  axial  dis¬ 
tance  to  be  considered. 

The  main  program  calls  upon  SPECIE  to  calculate  mole  fractions, 
molecular  weights,  and  mass  fluxes  at  all  radial  locations,  and 
calls  upon  THERMO  to  calculate  static  enthalpies  of  the  species. 

Total  enthalpies  are  calculated.  The  TRANSP  subprogram  is  called 
upon  to  calculate  a  turbulent  viscosity,  and  INTEG  to  calculate  ra¬ 
dius  as  a  function  of  stream  function,  and  also  to  calculate  certain 
coefficients  (AP,  APR)  needed  for  subsequent  finite  difference  cal¬ 
culations. 

The  main  program  writes  the  data  on  initial  profiles  as  a  function 
of  stream  function.  It  then  proceeds  in  a  stepwise  fashion  to  cal¬ 
culate  the  flow  variables  as  a  function  of  axial  location. 

After  the  writing  of  the  initial  data,  there  is  very  little  necessity 
or  possibility  for  operator  or  programmer  intervention  in  the  main 
program.  The  number  of  radial  locations,  NNN,  is  a  major  variable 
which  is  controlled  initially  and  internally.  The  variables  IK  and 
IA  define  the  locations  of  the  inner  and  outer  bounds  of  the  viscous 
flow  region.  In  the  present  application,  the  inner  bound  is  taken 
to  be  the  center-line.  The  outer  bound  is  always  IA  =  NNN-1 . 

Species  calculations  run  from  1  to  ISP,  which  is  defined  previously. 

A  variable  KUSK  defines  the  number  of iterationswhich  have  been  per¬ 
formed  at  an  axial  location.  The  iterations  stop  after  KUSK  =  3. 
Prior  to  the  start  of  calculations,  the  axial  position,  Z,  is  incre¬ 
mented.  IA  is  checked  to  see  if  the  maximum  number  of  streamlines 
permitted  has  been  reached.  If  IA  =  NILS  -  1,  a  doubling  of  stream 
function  step  size  is  incurred,  together  with  a  reduction  of  the 
number  of  steps  considered.  This  is  done  through  the  subroutine 
CHANGE.  In  the  present  program,  NILS  =  51  is  defined  as  data. 

NILS  must  be  odd.  The  resulting  profiles  are  written  after  CHANGE. 

Whether  or  not  a  change  occurs,  the  boundary  values  at  NNN  are  cal¬ 
culated  at  the  incremented  axial  position.  The  boundary  value  for 
velocity,  UEO  =  U(NNN)  never  changes  since  the  inviscid  flow  is  con¬ 
sidered  to  occur  at  constant  pressure.  The  velocities  at  the  other 
radial  locations  are  calculated  through  SETUP  and  GAUSS.  If  |(U(IA)  - 
UEO ) /UEO | >  TOL,  an  indicator,  IAP  is  set  up  to  tell  the  program  to 
increase  the  viscous  region  (this  is  not  done  until  after  all  iter¬ 
ations  are  complete).  The  species  concentrations  are  calculated  at 
NNN,  and  then  in  the  viscous  region  through  SETUP  and  GAUSS.  If 
the  relative  difference  in  molecular  weights  at  IA  and  NNN  exceeds 
TOL,  the  indicator  IAP  is  set.  The  total  enthalpy  at  NNN,  HT(NNN) 

=  HEO  is  constant  in  the  inviscid  flow.  In  the  viscous  region,  new 
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total  enthalpies,  HT(I)  are  calculated  through  SETUP  and  GAUSS. 

If  the  relative  difference  in  total  enthalpies  at  IA  and  NNN  exceeds 
TOL,  the  indicator  IAP  is  set. 

In  the  first  interation,  species  production  terms  at  the  previous  lo¬ 
cation  are  used.  In  the  second  and  third  iterations  the  species 
production  terms  at  the  midpoint  of  the  interval  are  approximated  by 
a  Taylors  series  expansion.  The  same  process  of  calculating  bound¬ 
ary  values,  then  viscous  flow, values,  of  the  flow  variables  is  fol¬ 
lowed.  Ini  between  iterations,  TRANSP  and  INTEG  are  called  upon  to 
obtain  improved  values  o?  turbulent  viscosity  and  the  terms  AP,  APR 
used  in  SETUP,  SETYP,  SETHP. 

After  the  third  iteration,  IAP  is  checked  to  see  if  the  flow  field 
should  be  expanded.  If  it  should  be  it  is  expanded.  If  it  should 
not,  previous  values  of  the  flow  variables  are  identified,  and  the 
program  returns  to  increment  another  axial  step.  The  results  are 
printed  every  50  steps,  as  defined  by  the  variable  KYS.  KYS  is  set 
to  50  near  the  start  of  the  main  program.  If  this  causes  too  much 
printing,  KYS  could  be  set  to  a  larger  value. 


Inputs  to  MAIN 


There  are  no  inputs  read  into  MAIN  by  tape,  cards,  etc. 


Expected  Outputs 


After  INIT,  the  MAI  !  program  writes  initial  profiles  of  selected 
flow  variables  as  a  function  of  stream  function. 

After  CHANGE,  the  MAIN  program  writes  resultant  profiles  of  selected 
flow  variables. 

Other  writing  is  done  by  other  subroutines. 
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2.4.2  SPECIE  (KK1,  KK2,  NU) 


Purpose  and  Uses  of  SPECIE 

SPECIE  is  capable  of  calculating  C1I(I)  (average  molecular 
weight),  X(J, I)  (mole  fractions  of  ISP  species)  and  ROU(I) 
(mass  flux)  ,  from  axial  location  KK1  to  KK2 

If  NU  «  0,  CM(I)  and  X(J,I)  are  transmitted 
If  NU  <0,  ROU(I)  is  also  transmitted 
If  NU  =*  0,  only  ROU(I)  is  transmitted 

Inputs 

There  are  no  system  inputs  on  cards  or  tape  to  SPECIE 
Expected  output 

There  are  no  system  outputs  on  cards  or  tape  from  SPECIE 
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2.4.3  THERMO  (KRIS,  JA,  IT) 


Summary  of  Operations  Performed  by  THERMO 


In  the  Initial  phases  of  the  progrrm,  JA  is  set  to  be  greater 
than  zero,  which  causes  THERMO  to  read  in  coefficients  for  calcu¬ 
lating  species  enthalpies.  When  JA  is  less  than  or  equal  to  zero, 
and  KRIS  equals  100,  species  enthalpies  are  calculated  at  T(IT) 
(temperature  allocation  IT);  if  KRIS  is  different  frcm  100,  tempera¬ 
ture  is  determined  at  a  static  enthalpy  HST(IT),  by  means  of  a  New- 
ton-Raphscn  process. 


System  Inputs 


A  total  of  ISP  sets  of  coefficients  A,  followed  by  ISP  coefficients 
B.  These  are,  respectively,  the  3OOrlO00°K  and  1000-5000'*K  coeffi¬ 
cients  ax  -  aQ  from  Table  V  of  NASA  SP-3001,  for  the  appropriate 
compounds.  For  fuel,  a  value  of  -23900R  calories  per  gram  mole  of 
CiQ  Hg2  has  been  estimated  to  be  the  enthalpy  of  fuel  at  477“F 
(estimated  boiling  point)  from  the  elements  at  298°K.  The  program 
also  reads  SENS,  heats  of  formation  at  298UK.  These  are  not  used 
in  the  present  program. 


System  Outputs 


The  coefficients  B  are  printed  after  being  read  in. 
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2.4.4  INIT 


Summary  of  Operations  Performed  by  INIT 

INIT  is  called  only  once.  It  reads  in  NNN,  the  number  of 
streamlines  considered,  then  NNN  cards  containing  UE(I) 
(velocity)  SE(I)  (temperature)  and  R(I)  (radius)  with  R(I) 
starting  at  zero  and  increasing  monotonically .  Initial 
profiles  of  species  mass  fractions  (YE(J,I))  are  then  read 
in  at  R(I) .  Stream  functions  are  calculated  for  the  initial 
profile,  and  the  data  are  interpolated  to  equal  intervals  of 
stream  function  with  the  interval  DPSI  =  PSI/(NNN-1),  where 
PSI(I)  is  the  value  of  stream  function  at  location  I. 

System  Inputs 

NNN  is  an  Integer  less  than  or  equal  to  NILS  (main  program)  . 
UE(I)  is  in  feet  per  second,  SE(E)  in  degrees  ,  Rankine, 

R(I)  in  inches.  Species  mass  fractions  are  introduced  per 
card  up  to  ISP,  for  each  R(I) 

System  Outputs 

Nothing  is  written  on  card  or  tape  from  INIT. 
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2.4.5  TRANS P  (NLOW,  NHIGH) 


Summary  of  Operations  Performed  by  TRANSP 

TRANSP  calculates  a  turbulent  viscosity,  UM,  according  to  Schetz' 
displacement  thickness  model,  as  described  in  the  report.  A  DATA 
statement  defines  RJ,  an  inner  jet  radius,  in  feet. 

System  Inputs  and  Outputs 

Nothing  is  read  or  written  by  TRANSP. 


2.4.6  INTEG  (NLOW,  NHIGH,  IN) 


Summary  of  Operations  Performed  by  INTEG 


INTEG  calculates  radial  values  (R1I)  as  a  function  of  stream  func¬ 
tion.  It  also  calculates  AP(I)  and  APR(I)  at  locations  NLOW  to 
NHIGH. 


If  IN  =  0,  AP  =  APR  =  SP 

If  IN  is  not  zero,  AP  =  .5  (SP  +  APR) 

where  SP  =  UM  *  R(I)  **2  *  ROU(I)/PSI (I) 

is  a  term  which  arises  in  the  transformation  of  equations  to 
stream  function  coordinates. 

System  Inputs  and  Outputs. 

No  information  is  read  to  or  written  from  INTEG. 
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2.4.7  SETUP  (RUP) 


Summary  of  Operations  Performed  by  SETUP 


SETUP  calculates  coefficients  for  solution  of  the  finite  difference 
equations  for  velocity.  These  are  transmitted  to  labelled  COMMON 
UPSET  and  used  by  GAUSS.  The  indicator  KOR  is  set  equal  to  -10. 


System  Inputs  and  Outputs 


There  are  no  INPUT/OUTPUT  Statements  in  SETUP. 
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2.4.8  SETYP  (EUP) 


Summary  of  Operations  Performed  by  SETYP 


SETYP  calculates  coefficients  for  solution  of  the  finite  difference 
equations  for  species  mass  fractions,  to  be  used  by  GAUSS.  The 
indicator  KOR  is  set  to  0.  A  DATA  statement  defines  CN  which  is  the 
inverse  of  the  turbulent  Prandtl  number.  CN  is  taken  to  be  1  in  the 
present  study. 


System  Inputs  and  Outputs 


There  are  no  INPUT/OUTPUT  Statements  in  SETYP. 
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2.4.9  SETHP  (RUP) 


Sunmary  of  Operations  Performed  by  SETHP 


SETHP  calculates  coefficients  for  solution  of  the 
equations  for  total  enthalpy,  to  be  used  by  GAUSS 
KOR  is  set  to  +  10.  A  DATA  statement  sets  CN  =  1 
the  present  program. 


System  Inputs  and  Outputs 


There  are  no  INPUT/OUTPUT  statements  in  SETHP. 


finite  difference 
The  indicator 
as  in  SETYP,  in 


2.4.10  GAUSS 


Summary  of  Operations  Performed  by  GAUSS 


Depending  on  whether  KOR  (see  previous  three  subroutines)  is  nega¬ 
tive,  zero,  or  positive,  GAUSS  solves  for  velocity  (U),  species 
mass  fractions  (Y(J,I))  or  total  enthalpy  (HT)  through  the  viscous 
flow  field.  When  KOR  =  I),  the  mole  fractions  (X(J,I))  and  average 
molecular  weights  (CM)  a\-e  also  determined.  When  KOR  is  positive, 
static  enthalpy  (HST)  ie  also  determined.  In  addition,  GAUSS  de¬ 
termines  the  outermost  location  at  which  the  flow  variables  U,  CM, 

HT  deviate  significantly  from  the  inviscid  external  flow.  This  is 
defined  to  be  IAP,  as  d.  cussed  in  MAIN. 

The  highest  location  at  which  significant  deviation  occurs  in  each 
branch  is  defined  to  be  NO(I),  I  =  1,  2,  3.  Quantities  NE(I),  de¬ 
fining  the  innermost  1<  nation  at  which  the  flow  variables  deviate 
significantly  from  an  j;  viscid  center  flow,  are  also  defined  but 
are  not  used  in  the  pr<  sent  program.  An  inner,  reacting,  inviscid 
flow  could  not  be  incl  i  led  in  the  program  without  some  modification 
(essentially  everythin/  done  in  the  stepwise  calculations  at  N^N 
would  have  to  be  done  it  the  inner  boundary).  An  inner  equilibrium 
flow  could  be  included  with  slight  modification. 

The  "significant  devi .  ion"  mentioned  above  is  defined  in  DATA  as  TOL, 
and  is  taken  to  be  1.  c  10“4  in  the  present  program. 


System  Inp  s  and  Outputs 
There  are  no  INPUT/O' '  PUT  statements  in  GAUSS. 
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2.4.11  CHANGE  (RUP) 


Summary  of  Operations  Performed  by  CHANGE 


CHANGE  is  called  upon  when  the  viscous  flow  field  has  expanded  to 
its  maximum  size.  The  value  of  DPSI  (stream  function  step  size) 
is  doubled,  and  the  number  of  intervals  is  halved.  In  the  present 
program,  NNN  is  reduced  from  51  to  26.  Values  of  the  flow  variables 
are  also  transferred  so  that  the  flow  variables  correspond  to  the 
proper  stream  function  values. 


System  Inputs  and  Outputs 


There  are  no  INPUT  Statements  in  CHANGE.  A  printed  output  indicates 
that  the  number  of  points  has  been  halved,  and  gives  values  of  cer¬ 
tain  variables  which  are  the  distance  downstream  from  the  start  in 
inches,  and  the  distance  interval  in  inches.  The  values  are  dup¬ 
licated  because  at  one  time,  DELX  was  changed  when  DPSI  was  changed. 
This  is  no  longer  done. 
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2.4.12  PRINT  (KUSK,  RUN) 


Summary  of  Operations  Performed  by  PRINT 

PRINT  lists  pertinent  variables  when  called  upon  by  the 
MAIN  program. 

System  Output 

The  first  data  listed  are 

Z(feet),  ZX(inches),  DELX(feet;.  DPSI,  NEL,  NEU,  IA, 
KUSK,  KUN 

Then  at  each  of  NNN  locations,  I  (radial  location),  PSI,  R, 
S (Temperature  in  °R)  AP,  APR,  HT  in  one  line,  followed  by 
Y(I , J)  (2  lines)  W(I,J)  (2  lines)  E(I,J)  (2  lines)  and  a 
final  line  with  YTOT,  CM,  ROU.  This  output  represents  the 
fundamental  output  for  the  program. 
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2.4.13  CHEMIE  (KRYS,  KRAS,  IT) 


Purpose  and  Uses  of  CHEMIE 

This  subroutine  calculates  quantities  necessary  to  calculate 
the  species  production  terms  for  inclusion  in  the  species 
conservation  equations.  It  also  reads  in  the  parameters 
for  calculating  reaction  ^ates.  If  (KRYS  +  KRAS)  is  not 
less  than  zero,  the  parameters  for  the  reactions  rates  are 
read  in  and  listed.  If  (KRYS  +  KRAS)  is  less  than  zero,  the 
reaction  rates  „  species  production,  and  linearization  terms 
are  calculated.  If  KRAS  is  greater  than  zero,  species  pro¬ 
duction  terms  (W)  are  calculated. 

If  KRAS  is  equal  to  zero,  linearized  terms  are  calculated. 
The  reactions  considered  are: 


H 

4* 

eg 

o 

- 

0 

+  OH 

(1) 

«2 

+ 

0 

= 

OH 

+  H 

(2) 

«2 

+ 

OH 

= 

H 

+  h2o 

(3) 

OH 

+ 

OH 

= 

h2° 

+  0 

(4) 

H2 

+ 

•  M 

= 

H  + 

H  +  M 

(5) 

H2° 

+ 

M 

= 

H  + 

OH  +  M 

(6) 

*2 

+ 

M 

= 

N  + 

N  +  M 

(7) 

°P 

+ 

M 

= 

0  + 

0  +  M 

(8) 

NO 

4- 

M 

= 

N  + 

0  +  M 

(9) 

n20 

+ 

M 

= 

N2 

f  0  4*  M 

(10) 

no2 

+ 

M 

= 

NO  +  0  +  M 

(11) 

CO 

+ 

OH 

= 

C02 

+  H 

(12) 

NO 

+ 

0 

N  + 

°2 

(13) 

NO 

+ 

N 

= 

N2  • 

t  0 

(14) 

NO 

+ 

°2 

= 

N°2 

+  0 

(15) 

NO 

+ 

NO 

= 

n2° 

+  0 

(16) 

C  H, 
n  : 

2n+2  + 

l°2 

- 

nCO  +  (n+1)  H2 

(17) 

System  Inputs  and  Outputs 


When  called  upon  near  the  start  of  the  program,  CHEMIE  reads 
data  for  calculating  chemical  reaction  rates.  These  are,  In 
order, 


CW  (molecular  weights) 

AU  (collision  terms  in  forward  rate  constants) 

CU  (collision  terms  in  backward  reaction  rate 

constants) 

E  (forward  rate  activation  energy/R) 

D  (backward  rate  activation  energy/R) 

GNUF(I,J)  (stoichiometric  coefficients  for  species 
I  in  forward  reaction  J) 

GNUB(I,J)  (stoichiometric  coefficients  for  species  I 
in  backward  reaction  J) 


AF,  BF, 
BB 

EN 


AB,  (respectively;  temperature  coefficient 


density  coefficient  for  forward  reaction 
similarly  for  reverse) 


(number  of  carbon  atoms  in  fuel  C 

n 


H2n+2) 
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Certain  of  the  above  are  listed  after  being  read  ir.,  namely, 
CW,  AU,  CU,  E,  D,  GNUF,  GNUB. 
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2.5  Usage  Instructions 
2.5.1  Preparation  of  Inputs 

Titles  and  Description  of  Inputs 

The  basic  inputs  to  the  program,  initial  profiles,  are  read 
by  the  subroutine  INIT.  The  titles  of  the  inputs,  read  in 
INIT,  are 

NNN  -  Number  of  inputs,  number  of  streamlines  considered 
initially. 

UE,  SE,  R  -  One  set  per  card. 

UE  is  specified  initial  velocity,  (ft/sec);  specified 
initial  temperature  (degrees  Rankin);  at  R,  radial  posit jon 
in  inches. 

YE(J,I)  A  total  of  ISP  species  mass  fractions,  arranged 
in  the  order,  H,  0,  H20,  OH,  0£,  H2,  N  ,  N,  NO,  N£0,  NO  ,  CO, 
C02,  fuel  at  R.  This  requires  two  cards  for  the  present  value 
of  ISP  =  14. 

The  data  are  read  on  canis . 

The  cards  must  be  read  in  the  sequency; 

NNN  Card 

UE,  SE,  R  -  Total  of  NNN  cards 

YE  -  One  set  of  cards  for  each  R  (2  NNN  cards  in  the 
present  program) . 

Limitations 

The  first  R  is  taken  to  be  zero.  R  must  increase  from 
card  to  card.  Negative  values  for  any  variable  would  have  no 
physical  meaning. 

Formats 


Rf AD (5, 100)  NNN 

READ15f  IQlMUEtUfSeq)  *R(U  »l«lfNNN) 
ftEA0(5,102H(YE(JtI)»J»l»ISP>»I«lfNNN) 


102  FORMAT! 7010.0) 


2  ( 


Relationship  of  Inputs  to  Outputs 

These  are  the  fundamental  starting  values  for  the  differ- 
entail  equations.  All  output  profiles  are  related  to  these. 

In  particular,  selected  flow  variables  are  printed  in  MAIN 
directly  after  INIT.  These  are  dependent  on  PSI ,  rather  than 
R,  as  are  all  subsequent  outputs  of  flow  variables.  Only  values 
at  the  end  points  should  be  the  same . 

2,5.2  Results  of  Operation 

Description  of  Results 

The  major  results  are  profiles  of  velocity,  temperature, 
and  species  mass  fraction  at  specified  axial  locations.  The 
output  comes  on  the  printer. 

Format  and  Content 

The  output  instruction  and  format  are  as  follows 

_ Z-ZX/12. _ 

WRITE  (6,30)  Z,ZX,DELX,DPSI  ,NEL,NEU, IKS,KUSK,KUN 

00  20  I -IK, II 

YTOT-O.O 

- no  io  - 

10  YT0T-VT0T+Y ( J, I > 

WRITE  (6,40)  I,PSI ( I  ) ,  H  (  I  )  ,  S  (  I  )  ,  U  ( I ) ,  AP  (  1 ),APR<  I ) ,HT(  I  ) 
WRITE  (6,50)  I,(Y( J,I >,J-l,tSP) 

WRITE  (6*50)  I,(W(J,D,J-1,ISP) 

WRITE  (6,50)  I , (EC J,I  ),j«l ,  ISP) 

WRITE  (6,60)  I , YT0T,CM( I ) , R0U( I  ) 

20 _ CONTINUE _ 

_ _  WR I  TE  (  6 ,60 )  KQ,UM _ 

RETURN 

C 

30  FORMAT  ( • 1 • , //✓ / , 4 ( 2X ,Fl 5 . 1 0 ) , 5 ( 3X , 1 5 ) , 7H  OUTPUT,///) 

40  FORMAT  ( / ,  2X , 14 , 7 ( 2X , 614. 7 ) ) 

30  FORMAT  ( 2X ,14,7 ( 2X , E14. 7 ) / , 6X ,7  t  2X , E14.  T ) ) 

60  FORMAT  (2X,  14,3 (2X,E 14.7))  " 


Use  of  Outputs 

The  outputs  contain  information  on  the  pollutant  concentra¬ 
tions,  and  other  concentrations.  Hence  the  use  of  output,  at 
the  end  of  the  afterburner,  will  be  to  predict  pollutant  con¬ 
centrations.  Intermediate  results  indicate  the  evolution  of  the 
pollutants . 
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2.6  Operating  Instructions 

2.6.1  Operating  Procedures 

Initiation  of  Program 

The  program  is  on  cards,  together  with  data.  For  a  par¬ 
ticular  computer,  proper  control  cards  must  be  introduced. 

Input  is  on  cards,  output  on  printer. 

Maintaining  Computer  Program  Operation 

The  program  will  continue  to  operate,  as  long  as  too  large 
a  step  size  is  not  used.  If  too  large  a  step  size  is  used,  the 
procedure  is  unstable,  and  is  characterized  by  logs  of  negative 
numbers,  due  to  temperature  going  negative.  The  program  must 
be  started  with  a  smaller  step  size  (DELX) . 

Termination  and  Restart 

Termination  is  achieved  when  Z  >  ZMAX.  If  interruption 
of  the  program  is  necessary,  some  subroutine  would  have  to  be 
written  to  identify  the  necessity  for  interruption  at  statement 
40,  MAIN,  store  program  and  all  data,  and  restart  at  the  same 
point . 

2.6.2  Operator  Inputs 

Initial  imputs,  other  than  the  basic  initial  profile  men¬ 
tioned  previously,  include  the  thermodynamic  and  chemical  data 
initially  read  in  THERMO  and  CHEMIE. 

The  input  statements  and  formats,  and  data  used  in  the 
present  version,  are  shown  below. 

In  THERMO 

- nro - READ  (5,7&UJ  I  “  - “ - 

- DO  19U” J*l  ♦  ISP  - - 

- RFAD  T5I  ZEUJ  TBTK7TJ )  f  K»1  *  6 ) - 

- WRITE  TBT7HPT  TUT  *  ;  J  )  i  K  « I  ,b ) - 

T9D - CONTINUE - 

-  - RF/nr_p5  ^770  J1 » J  a  17  T3P1 - 

GTTTmTn 

-  2U0 - WRrTE_T6T710T"N*"RTNJ  »"PSI  rWTf  HSTf  M)  *HTS»FA'»FBt  TP*T(  W) - 

- Off  TO  90  ~  '  ”  '  “ 

- c - 


FURMAT  T777»2K » I'4 ,  R  ( ZX  ,  E 1 3.3 )  *TX*4HTHER»  77TT 


- 270 FORMAT  T7TJT07TTJ - - 

~~  2RU - FORMAT'!  777771  IX»6EZU.H  7777777 
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The  Data  used  in  the  present  version  are  as  follows 


®*25O0OOM1  QL 


0.907775476-12  0.2913T19( 


luu 


-0.52546531E* 

3.371899466 

0.270821806’ 

0.284608496 

-0.33093421E- 

6.369161486 


13  0.358527676 
01-0.251672666* 
11-0*105767066 
01  0.419321166* 
11-0.967253726 
Ol-O. 13332532^* 
13-0.106283366 
or-.TT2437910f* 


•02  0.8 583 7 35 36—05-0* 829987 16E-08 


04 _ 

02-0.961193326-05  0.95122662E-08 
03 


•02  0.265031006— 05-0 c 97688 341 E-0 9 
04 

•03 .29647566006-06-.  3246404906-09 


.22?095l2006- 11.97647894006 


.34344563006 

•3721252300E* 

TTfijTmzWF 

>772 1684100E* 

♦2l701000Qgr 

-.1628070106- 


01.22234297606 
-ll. 28647685006 


02 . 9692246700E— 05-.  786336390E-08 


04 


01-.  1 1 1676340E— 04.6958316500E-08 


04 


•02.6714B97500E— 05-.974277190E-08 


04 


12-. 1436350806 
01.1037811500 
11-. 4835260206 

^”,"23*9000  Q0QE 


-•  107339380E— 04.6345917500E-08 


i 


Blank  Card 


0.25372567E  Ol-O 
-0 . 5 5743608 E- 15  0 
0. 26707532E  0l  0 
-0.61973  5'6'BE- 1 4-(5" 
0.288955446  01  10 


»25470497E 
,18422l90t- 
>292300076 
►  303171156* 


nitroirf 


05 


•04-0 .88 
05 

•02-0, 8535 15706-06  0. 1 1790853E-09 


08  0.596436216-11 


r38811792 


-0.  334602045-14-0.  tl9Z79TB 


.24422261006  01. 1227618700E-03-. 8499271906-07, 2140083000E-10 


.T152936O0OE  01.1405995500 
-,7372078306-14.98522040006  04 


06, 1062820900E-09 


-.1584970106-10.1 
,46 1 3921 900 E  Ol.i 


1094854TOF— 05. 2081842500E-09 


■  01 .1 55255670? 

a- STTfi 


Lil  V !•>  I .T.l •: »1 i)  PM SjMtfMiM »l.SS’Z5t%fTV5»Tdgil- 


-.1674298606-13-7 


Rlank  Card 


_ -v.23  90nnom)E  Q5  _ _ 

59556 . f  -57T97. 9  9312 . 3  ~ 

716007  19490.  7007.5  ^7 


READ ( 5.606 )  EN 


Data  are 


“  1.008  TS75T? 
14,00b  W.008  ~  44.016 
.720001  I5.33000E  IT. 620001 


28.016  CW 
142.286  CW 
A 


The  GNUB  are 


8 
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SECTION  III 

LISTING  OF  PROGRAM 


A  listing  of  the  program  and  data  input  are  presented  in 
this  section 
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